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ABSTRACT 

In this paper, we analyze asymptotic delay-throughput trade-offs 
in mobile ad-hoc networks comprising heterogeneous nodes with 
restricted mobility. We show that node spatial heterogeneity has 
the ability to drastically improve upon existing scaling laws estab- 
lished under the assumption that nodes are identical and uniformly 
visit the entire network area. In particular, we consider the situation 
in which each node moves around its own home-point according to 
a restricted mobility process which results into a spatial stationary 
distribution that decays as a power law of exponent S with the dis- 
tance from the home-point. For such restricted mobility model, we 
propose a novel class of scheduling and routing schemes, which 
significantly outperforms all delay-throughput results previously 
obtained in the case of identical nodes. In particular, for 5 = 2 
it is possible to achieve almost constant delay and almost constant 
per-node throughput (except for a poly-logarithmic factor) as the 
number of nodes increases, even without resorting to sophisticated 
coding or signal processing techniques. 

1. INTRODUCTION 

Over the last decade we have seen a flurry of theoretical studies 
aimed at establishing fundamental scaling laws of ad-hoc networks 
as the number of nodes increases. Gupta and Kumar first consid- 
ered the case of n static nodes and n random source-destination 
(S-D) pairs, obtaining the disheartening result that the maximum 
per-node throughput decays at least as l/\/n [1]. 

In contrast to static networks, Gossglauser and Tse [2] have 
shown that a constant per-node throughput can be achieved in mo- 
bile ad-hoc networks by exploiting the store-cany-forward com- 
munication paradigm, i.e., by allowing nodes to store the data and 
physically carry them while moving around the network. The result 
in [2] was proven under the assumption that nodes independently 
move according to a generic, ergodic mobility process which re- 
sults, for each node, into a uniform stationary distribution over the 
space. This mobility model is actually a generous one, as it allows 
each node to equally come in contact with any other node, achiev- 
ing a full, homogeneous mixing. 

In practical cases, however, the mobility pattern of individual 
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nodes is expected to be restricted over the network area, as users 
spend most of the time in proximity of a few preferred places [3], 
and rarely go outside their region of habit. This observation has al- 
ready motivated some researchers to study the impact of restricted 
mobility models. In [4] a one-dimensional mobility model is con- 
sidered, in which each node uniformly visits a randomly chosen 
great circle on the unit sphere, obtaining again a constant through- 
put. In [5,6] the authors consider a two-dimensional, restricted mo- 
bility model which produces, for each node, a rotationally invariant 
spatial distribution centered at a home-point uniformly chosen in 
the area; the resulting throughput varies with continuity in between 
the two extreme cases of static nodes (Gupta-Kumar) and fully 
mobile nodes (Grossglauser-Tse), depending on how the physi- 
cal network extension scales with respect to the average distance 
reached by the nodes from their home-point. This result confirms 
that throughput is maximized when the nodes span the entire ex- 
tension of the network area. However, the authors of [5,6] have not 
analyzed the delay under their restricted mobility model. 

Driven by the optimality of the homogeneous mixing assump- 
tion in terms of throughput, many authors have analyzed asymp- 
totic delay-throughput trade-offs under the same assumption. This 
choice is also motivated by the fact that the most popular mobil- 
ity models adopted in the literature (such as random walk, random 
way-point) produce a uniform stationary distribution over the area[j 
Indeed, when considering also the data transfer delay, the precise 
details on how the nodes move become important. Several mobil- 
ity models have been analyzed, ranging from the simple reshuffling 
model [7], to the Brownian motion [8], and variants of random walk 
and random way-point [9, 10]. In all of these studies, nodes have 
been assumed to be identical and fully mobile, i.e., their trajecto- 
ries 'fill the space' over time, uniformly visiting the entire network 
area. 

Starting from a different perspective, essentially aimed at re- 
ducing the large delays encountered in sparse, intermittently con- 
nected mobile ad hoc networks, many researchers have already 
proposed distributed routing protocols exploiting the novel store- 
carry-forward communication paradigm, within the context of 
Delay Tolerant Networking (DTN). Some of them have already 
pointed out that node heterogeneity in terms of spatial locality or 
inter-contact times can be very beneficial to improve end-to-end de- 
lays. In particular, the history of past encounters [11,12] or the ex- 
plicit dissemination of information about the mobility pattern of the 
nodes [18, 19] can significantly improve the message delivery de- 
lay without resorting to flooding-based approaches which are very 
wasteful of resources. 



the stationary distribution of a node under the random way-point 
is uniform in the absence of border effects, such as on the surface 
of a sphere or a torus. 



However, the feasible performance gains, in terms of both 
throughput and delay, that can be achieved by exploiting node spa- 
tial heterogeneity, as well as the resulting scaling laws in a network 
with increasing number of nodes have not been investigated so far. 
In this paper, we bridge the theoretical analysis of fundamental 
scaling laws of mobile ad-hoc networks with the insights already 
gained through practical protocol development. By so doing, we 
provide a theoretical foundation to the design of intelligent routing 
schemes which exploit the spatial heterogeneity of nodes, analyt- 
ically showing the potential of such schemes in terms of delay- 
throughput trade-offs. 

In particular, we consider a restricted mobility model similar 
to the one introduced in [5], in which each node moves around a 
home-point randomly chosen in the area. Nodes move indepen- 
dently of each other, but they are not identical, because each node 
is characterized by a different home-point. We consider an ergodic 
mobility process which produces, for each node, a spatial stationary 
distribution which is rotationally invariant around the home point, 
and decays as d^*, where d is the distance from the home-point, 
and (5 > (the distribution is properly normalized so that its in- 
tegral over the area is one). In order that each node achieves the 
above spatial distribution we have considered, for simplicity, the 
reshuffling model, usually referred to as i.i.d. mobility model. 

The family of mobile networks that we consider comprises, as 
a special case, the Grossglauser-Tse scenario in which nodes are 
fully mobile (5 = 0) and therefore indistinguishable from each 
other, as well as the limiting case of static nodes considered by 
Gupta-Kumar (S — > cxa). We identify a class of scheduling-routing 
schemes whose performance in terms of throughput and delay ex- 
hibits an intriguing behavior as we vary S (see Figure [T]!. In par- 
ticular, for S — 2 our scheme achieves near-optimal results, i.e., 
almost constant throughput and almost constant delay (except for a 
poly-logarithmic factor), and, over a wide range of values for 5, it 
significantly improves over existing bounds derived under the as- 
sumption of identical nodes. 

We emphasize that our schemes do not exploit any sophisticated 
technique like the ones that have recently been proposed to improve 
delay-throughput trade-offs in the basic cases of Gupta-Kumar and 
Grossglauser-Tse, such as hierarchical cooperation with MIMO 
communications [13] and source coding [14], respectively. The 
purpose of our work is not to establish optimal information theo- 
retic results, but to show that there is an additional dimension to 
be exploited, i.e., node spatial heterogeneity, which has been so far 
neglected by theoretical studies aimed at establishing fundamen- 
tal scaling laws of mobile ad-hoc networks. For this reason we 
will maintain the basic system assumptions originally introduced 
by Gupta-Kumar. More sophisticated techniques can be added to 
our scheme, and can further improve the bounds presented here. 

2. SYSTEM ASSUMPTIONS 

2.1 Mobility model 

We consider a network composed of n nodes moving over a 
square region O of area n with wrap-around conditions (i.e., a 
torus), to avoid border effects. Note that, under this assumption, 
the node density over the area remains constant as we increase n, 
equal to 1 . 

Time is divided into slots of equal duration, normalized to 1. 
We consider a two-dimensional i.i.d. mobility model, according to 
which the positions of the nodes are totally reshuffled after each 
slot, independently from slot to slot and among the nodes. At the 
beginning of each slot, a node jumps in zero time to a new position, 
and remains in the new position for the entire duration of a slot. 



Although the i.i.d. mobility model may appear to be unrealistic, 
it has been widely adopted in the literature because of its mathemat- 
ical tractability. We have adopted it in our work especially because 
it allows to model in a straightforward way the situation in which 
the stationary distribution of a node over the space is not uniform, 
which is more difficult to obtain using other mobility models com- 
monly adopted in the literature. 

Similarly to previous work, we consider two time-scales of node 
mobility: 

Fast Mobility: the mobility of nodes takes place at the same time- 
scale as packet transmission. Therefore, when two nodes 
wish to communicate, they have only one time slot at their 
disposal, after which the two nodes separate from each other. 
As a consequence, only single-hop transmissions can occur. 

Slow Mobility: the mobility of nodes is sufficiently slow that a 
node can send a packet over multiple hops to reach another 
node. This situation is usually modelled by redefining the 
time slot as the duration of the 'coherence interval' during 
which the positions of the nodes can be considered to be 
static. The topology is reshuffled after each slojj, but sev- 
eral transmissions can occur during a slot because the packet 
transmission time can be set much smaller than the duration 
of a slot. 

In our work we will consider both fast and slow mobility, but we 
warn the reader that the results obtained under these two scenarios 
are not directly comparable, because the definition of time slot is 
different. 

Let Xi{t) denote the position of node i at time t (t is 
an integer denoting the slot sequence number) and X(t) = 
{Xi{t),X2{t) ■ ■ ■ Xn{t)) be the vector of nodes' positions; we de- 
fine by dij (t) the distance between mobile i and mobile j at time 
t,Le^d,j{t) = \\X,{t)-X,{t)\\ 

Each node i is characterized by a home-point Hi, which is uni- 
formly and independently selected over the area. The collection of 
the nodes home-points is denoted by vector H = {Hi , H2 ■ ■ ■ Hn) 
and does not change over time, although it can be different for each 
network instance. We define by d^ the distance between the home- 
points of nodes i and j, i.e., df^ = \\Hi — Hj\\. 

The spatial stationary distribution of a node is assumed to be ro- 
tationally invariant with respect to the home-point, and thus can be 
described by a generic, non increasing function (/>(d) of the dis- 
tance d from the the home-point. In this paper we assume that 
4>{d) decays as a power-law of exponent S, i.e., (j>{d) ~ d~ , 
with S > 0. This choice is supported by a number of measure- 
ments papers which have found power-laws to be quite ubiquitous 
in experimental traces related to both human and vehicular mobil- 
ity [16-19]. For example, in [17] authors analyze a large corporate 
wireless network and find that the fraction of time spent by users in 
association to a given access point exhibits a power law. 

To avoid divergence problems in proximity of the home-point, 
we take function s{d) — min(l, d~ ), and normalize it so as to 



In the slow mobility case, the reshuffling (i.i.d.) mobility model 
can be justified assuming that devices are disconnected from the 
network (e.g., switched off) while travelling from one place to an- 
other. 

^^ Given any two points Xi — {xi,yi)€^0 and 
X2 = {x2,y2)&0 we formally define their distance 
over the torus surface as: |jXi(t) — X2(i)|| = 

min„^„g{_^_o_^} ^(xi +u- X2Y + (yi + i; - 3/2)^ 



obtain a proper probability density function over the network area: 

s{d) 
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The value of the normalization constant G — JJ„ s{d) can be ap- 
proximated, in order sense Q by the following integral in polar co- 
ordinates: 



G = e 



if AO I min(l,p"'')pdp j 



We obtain that G is finite for any d > 2. For < 5 < 2 we 
have G — Q{n 2 ). For the special value = 2 we have G = 
0(log n). Note that 5 = Q leads to a uniform distribution over the 
space, whereas letting 5 go to infinity we obtain the same behavior 
as that of a static network. 

2.2 Communication model 

To account for interference among simultaneous transmissions, 
we adopt the protocol model introduced in [l]|j. Nodes employ 
a common range R for all transmissions which occur in the same 
time slot; equivalently, they employ a common power level, i.e., no 
power adaptation mechanism is used. A transmission from node i 
to node j using transmission range R can be successfully received 
at node j if and only if the following two conditions hold: 

1. the distance between i and j is smaller than or equal to R, 
i.e., dy(i) < R 

2. for every other node k simultaneously transmitting, 
dkj{t) > {1 + A) R, being A a guard factor. 

Transmissions occur at fixed rate which is normalized to 1. We 
assume that a single copy of each packet is present in the network 
at any time, i.e., data units are not broadcasted nor replicated, and 
nodes do not keep copies of previously received packets in their 
buffer. 

2.3 Traffic model 

Similarly to previous work we consider permutation traffic pat- 
terns in which n randomly selected source-destination pairs ex- 
change traffic at rate A. Source-destination pairs are selected is 
such a way that every node is origin and destination of a single traf- 
fic flow with average rate A. We further assume that, for each pair, 
the distance between the source home-point and the destination 
home-point is Q{^/n), i.e., we consider the worst case in which all 
connections are established among nodes having 'far away' home- 
points. 

2.4 Throughput and delay 

We use the following definitions of asymptotic throughput and 
delay. Let Li (T) be the number of packets delivered to the desti- 
nation of node i in the time interval [0, T]. The delay of a packet 
is the time it takes for the packet to reach the destination after it 

*Given two functions f{n) > and g{n) > 0: f(n)= o{g{n)) 
means lim„^oo f{n)/g{n) — 0; f{n) = 0(g{n)) means 
limsup„^ /(n)/g(n) = c < 00; f{n) = io{g{n)) is equiv- 
alent to g{n) — o{f{n)); f{n) — Q{g{n)) is equivalent to g{n) — 
Oifin)); f{n) = e(<?(n)) means f{n) = 0{g{n)) and 5(71) = 
0(f{n)); at last f{n) x g(n) means liin„^oo f{n)/g{n) = 1 
^The protocol model has been proven to be pessimistic with respect 
to the physical model employing power control (see Theorem 4.1, 
pag. 174 in [20]). Thus the results obtained in this paper can be 
regarded as lower bounds of the network performance achievable 
under the physical model employing power control. 



leaves the source. Let Di{t) be the sum of the delays experiences 
by all packets successfully delivered to the destination of node i 
in the time interval [0, T]. We say that an asymptotic throughput 
A and an asymptotic delay D per S-D pair are feasible if there is 
an no such that for any n > no there exists a scheduling/routing 
scheme for which both the following properties hold 
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Equivalently, we say in this case that the network sustains an ag- 
gregate throughput A = nX. We will also adopt the simple Power 
Function [21], defined as the ratio X/D, to characterize the system 
performance by a single metric. 

3. PREVIOUS RESULTS FOR THE I.I.D. 
MOBILITY MODEL 

To avoid confusion, we limit ourselves to reporting existing scal- 
ing laws obtained for the i.i.d. mobility model. We emphasize that 
all of the following results have been derived under the assumption 
that the nodes' spatial distribution is uniform over the area. 

In [15] the authors analyze throughput/delay trade-offs in the 
case of fast mobility, with or without packet replication. The fol- 
lowing general trade-off is established: 



X^Oi^ 

n 



(1) 



In particular the two-hop scheme of Grossglauser-Tse is proven 
to incur a delay D — 0(n) while guaranteeing a per-node through- 
put A = 0(1). To improve delay, two different schemes exploit- 
ing packet redundancy are proposed: the first is still based on 
two hops, and achieves D — O(v'n) and A = 0{l/^/n); the 
second employs multiple hops and achieves better delay perfor- 
mance D — O(logn) while sacrificing the per-node throughput 
A = 0(l/(nlogn)). 

The i.i.d. mobility model in the case of slow nodes has been 
studied in [7,22]. In [7] a class of scheduling-routing schemes that 

achieves A = O ( •y/D/nlog"'^" n) was devised. A better trade- 
off was obtained in [22], where the authors show that necessarily 
the following relation must hold: 




and they propose a scheme that approaches this bound up to a poly- 
logarithmic factor, even in the case of constant delay. 

Recently it has been shown in [14] that previous results can be 
further improved by encoding transmitted information at sources. 
A class of joint coding-scheduling-routing schemes is introduced 
which achieves a trade-off A — Q{\/ D/n) in the case of fast mo- 
biles, and A = B( yDjn) in the case of slow mobiles. 

4. SUMMARY OF RESULTS 

A graphical representation of our results is reported in Figure [T] 
and Figure[2l Here we discuss only the fast mobility case. Results 
for the slow mobility case are similar and will be presented in Sec- 
tion |7] In Figure [T] we have plotted, as a function of 5, the best 
Power XjD achievable by the class of scheduling-routing schemes 
introduced in this paper. We have employed a log^ scale in the ver- 
tical axis so as to show the asymptotic order in n. Note that in this 
scale we can neglect the impact of poly-logarithmic factors. 
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Figure 1: Power \/D (in log scale) as a function of 5 (fast mo- 
bility case). 



We observe that the maximum Power according to our scheme is 
achieved when 5 = 2: for this particular value both A and D remain 
almost constant with n, resulting in a Power which scales as n'^ = 1 
except for a poly-logarithmic factor Hence, the performance of our 
scheme for 5 = 2 is very close (in order sense) to the best result 
that one can think of. For 1 < 5 < 2 the optimal Power has the 

expression n <'-'' , whereas for 2 < (5 < 3 the Power is given 
by n^~ . For all other values of 5 the best Power is equal to —1. 

For 5 > 2 no delay-throughput trade-offs are possible employ- 
ing our schemes, since throughput is maximized jointly with the 
minimization of delay. 

A wide range of delay-throughput trade-offs is instead possible, 
within our class of scheduling-routing schemes, for 5 < 2. How- 
ever, for 1 < 5 < 2, such trade-offs may result in a sub-optimal 
Power This is illustrated in Figure|2l in which we have plotted (in 
logjj-log^ scale) all feasible combinations of A and D for various 
values of 5. For a given 5, the feasible combinations of A and D lie 
on a line departing from the common point A = 1, D = n. The 
little circles, which lie on the curve D — A^^, denote, for each 
considered 5, the smallest achievable delay (corresponding also to 
the smallest throughput), and correspond to the point in which the 
Power function is maximum (as reported in Figure[TJ . Notice that, 
for (5 < 1, all points lie on the curve D = n\, having common 
Power equal to —1. Notice that this result agrees with the general 
trade-off (TJ derived in [15] for fast mobiles uniformly distributed 
in the space. Restricted mobility starts to be beneficial in terms of 
delay-throughput trade-offs when 5 > 1. In particular, as 5 ap- 
proaches 2, it is possible to push down the delay towards 1 at the 
expense of smaller and smaller degradation of the throughput. For 
example, for 5 — 1.95 one can achieve a delay close to 1 with just 
a little penalty in throughput. We have also shown on the graph the 
curve D — n}? , which represents the best trade-off available so 
far for the fast mobility case, obtained in [14]. We observe that 
our schemes performs better as soon as 5 > 4/3, even without 
resorting to coding techniques. 

We now provide an intuitive explanation of the above results. 
Our schemes exploit node heterogeneity by an intelligent selection 
of relay nodes based on the location of their home-pointsQ Essen- 
tially, we adopt a geographical routing strategy combined with a 
divide-and-conquer technique. Data are forwarded along a chain 
of relay nodes whose home-points progressively close in on the 
home-point of the destination. At each step, the distance between 
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we assume that each node knows the location of its home-point 
and of the home-points of nodes falling in its transmission range. 



Figure 2: Delay-Throughput scaling trade-off (in log-log scale) 
for different values of 5 between and 2 (fast mobility case). 
The marks on the axes represent the orders asymptotically in 



the home-point of the next-hop relay and the home-point of the 
target is halved, guaranteeing that the message is delivered to the 
destination in log n steps. Now, at a given point in space the density 
of nodes whose home-points are at distance 0{Z) from the point- 
grows as Z^Z"^ jG. Therefore, for 5 < 2 a node most frequently 
gets in contact with nodes whose home-point reside far way in the 
network area. As a consequence, for 5 < 2 the critical step for 
the system performance is the last one, in which the packet has to 
advance by the minimum distance along the chain of home-points. 
Instead, for <5 > 2 a node typically gets in contact with nodes whose 
home-point are close to the home-point of the node. As a conse- 
quence, for 5 > 2 the critical step for the system performance is 
the first one, in which the packet has to advance by the maximum 
distance along the chain of home-points. The value (5 = 2 is the 
unique exponent at which the home-points of the nodes encoun- 
tered by a given node are equally distributed at all distance scales. 
As a consequence, no step is critical and the system achieves the 
optimal performance both in terms of throughput and delay. 

We remark that our system presents an interesting analogy with 
the problem of navigating small-world graphs using decentralized 
algorithms employing only local contact information. In particu- 
lar we mention here a well known result due to J. Kleinberg [23], 
who studied a 2-dimensional lattice enriched by random shortcuts 
according to a probability which decays as a power law with the 
distance between the connected vertices. The number of hops re- 
quired to reach an arbitrary destination exhibits a behavior similar 
to that in Figure [T] as a function of the power-law exponent. In 
particular, a unique value of the exponent (equal to 2) allows to 
navigate the graph in log n hops. 

5. FAST MOBILITY: SCHEDULING- 
ROUTING SCHEMES 

In this section we describe the scheduling-routing schemes that 
we have devised for our system, in the case of fast mobility. They 
form a family of schemes because one parameter, Zo, is free, and 



can be specified within a particular scheme so as to achieve a de- 
sired delay-tiiroughput trade-off. 

Before going on we premise a useful concentration result widely 
used in previous work. 

Lemma 1 . Let T be a tessellation of O, whose elements Ti 
have area \Ti\. The number of home-points falling within each Ti 
is, with high probability, O(jTij), uniformly over O, as long as 

lT,l = n(logn),Vi. 

We do not repeat the proof of this lemma, which is based on a 
standard application of the Chernoff bound (see [24]). 

5.1 Routing schemes 

As already mentioned, we propose a bisection technique that 
makes messages advance along a chain of relay nodes whose home 
points become progressively closer to the home point of the des- 
tination. We stop the bisection when the home-point of the target 
node falls within distance Zo from the home-point of the last relay 
node. 

To formalize this idea, suppose that node a handles a message 
m directed to d. Node a first computes d^ = \\Ha — Hd\\', if 
dad "£ Zo node a directly forwards the message to node d at the 
first transmission opportunity among the two nodes (as dictated by 
the scheduling scheme). 

If d^ > Zo, node a compares d^ with the set of thresh- 
olds {Zi = 2'Zo, i G N}, finding the i* > 1 such that, 
2' ~^Zo < d^d < 2' Zq. In this case node a is allowed to for- 
ward the message to any node b whose home-point satisfies the 
following geometric constraint with respect to the home-point of d: 
2^''^Zo < d^ < 7 2'*"^^o, where 1/2 < 7 < 1 is a constant 
introduced to simplify the analysis of our schemes. For the sake 
of concreteness, in the following we will consider the special case 
of 7 = 3/4. However, notice that asymptotic results are, in order 
sense, insensitive to 7. While being handled by a, we say that mes- 
sage m is in step i* of the routing algorithm. Note that steps are 
counted backward starting from the destination, and that message 
m has to go through i* relay nodes before being delivered to the 
destination. 

As an example, consider the case illustrated in Figure[3] in which 
source node a wants to send a message to node d. Distance df^ 
between the home-points of a and d satisfies 4Zo < df^ < 8Z0, 
hence we have i* — 3. This means that message m is in step 3 of 
the routing algorithm, thus it has to go through 3 relay nodes before 
being delivered to the destination. Since in this case 2' ^'^Zo = 
2Zo and |2' ~^Zo — SZq, the home-points of the nodes b to 
which a can forward the message have to lie in the ring 2Zo < 
d^ < 3Zo, which is depicted in Figure[3]as a shaded region. 

5.2 Scheduling schemes 

A scheduling scheme is in charge of selecting, at each slot, the 
set of (non-interfering) transmitter-receiver pairs to be enabled in 
the network, as well as the message to be transmitted over each 
enabled pair. In our family of schemes, the transmission range em- 
ployed by a transmitter depends on the routing step reached by the 
message to be sent. To better pack simultaneous transmissions, 
and thus maximize the network throughput, our scheduling policy 
selects, in a given slot, transmissions having homogeneous trans- 
mission ranges. This is done by selecting, at the beginning of each 
slot, a step i according to an assigned probability distribution pf . 
Then the slot is reserved only to messages which are in step i of 
the routing algorithm, i.e., to messages currently stored at nodes 
whose home-points are at distances ranging between 2^^^Zo and 
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Figure 3: Example of routing of a message from a to d. Dots 
denote home-points of the nodes. 
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Figure 4: Example of non-interfering subset of squarelets with 

M = 9. 



2' Zo from the home-point of the destination (for the particular case 
j = 0, distances are between and Zo). 

A common transmission range Ri is employed by all commu- 
nications occurring during a slot devoted to step i. Once step i 
has been selected, the domain O is divided into squarelets of area 
Ai = R^, forming a regular square tessellation. According to the 
protocol model, at most one transmission can be enabled in each 
squarelet. Moreover, one can easily construct AI = 0(1) sub- 
sets of regularly spaced, non-interfering squarelets (for example, 
M — 9 assuming a protocol model with A = 0, see Figure |4). 
Each subset can then be enabled in one out of M slots, guaran- 
teeing fairness among all squarelets and absence of interference 
among concurrent transmissions. 

Within a given squarelet 5, the specific transmitter-receiver pair 
(a, b) to be enabled is selected as follows. 

In slots devoted to the last step (i — 0), first all pairs (a, d) 
residing in 5 and such that: i) df^ < Zo; ii) a has some message m 
destined to d, are first classified as eligible for transmission. Then, 
if the set of eligible pairs is not empty, the squarelet 5 is declared 
as active and one pair is randomly selected for transmission. Notice 






(a) Examples of locations of Ha at 
routing step 3, given Hd 



(b) Examples of locations of Hd at routing step 
3, given Ha- The shaded regions around H'd and 
H'd denote the corresponding locations of Hb 



(c) The shaded region denotes the union of all 
possible locations of Hi, at routing step 3 



Figure 5: Geometric considerations supporting tlie proposed scheduling sctiemes 



that, provided that Zq — Q{^/\ogn), by Lemma [T] each node has 
w.h.p. G(Zo) destinations d. 

In slots devoted to step z > 0, first each node a residing in 5 
selects a message m from its buffer, which has reached step i (if 
there is any); let dm be the final destination of message m. All pairs 
of nodes (a, b) residing in 5 and such that: i) a has selected one 
message m; ii) b can act as (i — l)-th relay for message m, are first 
classified as eligible for transmission. Then, if the set of eligible 
pairs is not empty, the squarelet 5 is declared as active and one 
pair is randomly selected for transmission. We further assume that 
message m is selected by node a according to a FIFO scheduling 
policy. 

Thus, logically each node is equipped with a FIFO queue for 
each step i, i > 1, plus one queue per destination for step 0, in 
which messages are sent directly to the target node (see Figure[6}. 

An illustration of how our scheduling scheme works is reported 
in Fig. \5\ Suppose, as an example, that the current slot has been 
assigned to step i* = 3. Then we have a situation similar to the 
one already considered in Figure [3] while describing the routing 
scheme. Figure |5(a)| shows, for a given destination home-point 
Hd, two possible locations H'a and H" for the home-point of the 
transmitting node a. In particular, H'a (H'a) is characterized by the 
smallest (largest) distance from Hd- If we now take the point of 
view of a given transmitting node a, we obtain the symmetric situ- 
ation depicted in Figure [5(b)| in which H'd (H'd) are located at the 
smallest (largest) distance from Ha- The shaded ring around H'd 
(H'd) denotes the possible location of the home-point Hi, of feasi- 
ble relay nodes for the corresponding destination. Notice that, by 
selecting the head-of-the-line message, we constrain ourselves to a 
given destination home-point, and thus we can only use relay nodes 
whose home-points lie within a ring similar to the ones depicted in 

Figure|5(b)] 

In Figure [5(c)] we have shown the union of all of the possible 
rings of the type shown in Figure |5(b)[ as we let the location of 
Hd vary. The union region comprises all nodes b satisfying -^ < 
dab ^ ^Zi, for the generic step i > 1. 



Notice again that, provided that Zi = Q,{^/\ogn), by Lemma[T] 
there are w.h.p. Q{Zf) nodes b which may potentially act as relay 
node, and all of them are at distance OiZA from a. Moreover, the 



same is true if we constrain ourselves to a given destination home- 
point (the area of any shaded ring as in Figure [5(b)] is 0(Zf ), and 
all of its point are at distance 0{Zi) from a). This observation 
already suggests that our policy can obtain the same asymptotic 
performance (in order sense) as the one achievable by a more ag- 
gressive scheduling policy which selects the message to transmit 
based on the availability of next-step relays in the same squarelet 
of the transmitter. This intuition will be confirmed later in our anal- 
ysis. 

In the following, we will always assume that Zq (and thus all 
Zi's) are n{^\ogn), so that we can apply the concentration result 
in Lemma [T] 



6. ANALYSIS 
SCHEMES 



OF THE PROPOSED 



6.1 Design considerations 

Given that step i has been selected at the beginning of a slot, 
by construction the number Ni of parallel transmissions that can 
occur in the network during the slot equals the number of active 
squarelets. On average we have: 



E\m = 



MA, 



Pr{active squarelet | i} 



(2) 



where n is the total network area. A/ is a finite constant account- 
ing for interference (see |5.2t and Pr{active squarelet | i} is the 
probability that a generic squarelet S is active at step i- 

From the discussion in Section [5!2l a squarelet 5 is active if the 
following two conditions hold: i) at least one pair (a, b) of nodes 
is found within 5, such that b can act at as (i — l)-th relay node for 
a (i > 1), i.e., -^ < dab '^ ^^i (in the case i = 0, b must be 
a destination node satisfying -^ < dfl < f ^o); ii) node a has 
a message to transmit to node b (recall that, for i > 1, only the 
head-of-the-line message in the queue associated to step i can be 
transmitted). 

We observe that the probability that condition i) above holds, 
denoted by Pr{populated squarelet | i}, depends only on the ge- 
ometry of nodes and on the mobility process, not on the traffic (i.e.. 



queues dynamics). In general we have 



class of scheduling-routing schemes are the following: 



Pr{active squarelet | i} < Pr{populated squarelet | i} 

The choice of Ai is critical. The selection of a too large value for 
Ai leads to a suboptimal exploitation of spatial reuse (thus caus- 
ing throughput degradation), without being effective in reducing 
the delivery delay. This is because of the contention delay among 
many eligible transmitter-receiver pairs (a, 6) residing in squarelet 
S (recall that only one of them can be enabled), which offsets the 
advantage of reaching a more distant receiver in a single slot. 

On the other hand, the selection of a too small value for Ai is 
ineffective in terms of throughput (and also in terms of delay). This 
is because the potential benefit of increasing the spatial reuse is 
countered by the fact that the fraction of active squarelets tends to 
0. 

Therefore, a reasonable design choice is to minimize the 
squarelet size Ai under the constraint that: 

lim inf Pr{populated squarelet | i} > 

n — >oo 

Although this criterion provides only an upper bound to the prob- 
ability that a squarelet is active, we will later show in Section [63] 
that the network queues can be loaded, under our schemes, in such 
a way that the probability that node a indeed has a message to trans- 
mit to 6 is also w.h.p. greater than zero. 

We will further assume that \/]47 = Ri — 0{Zi); i.e., the trans- 
mission range used at step i should be smaller than or equal (in 
order sense) to the distance between the home-points of the trans- 
mitter and the receiver. This because it does not make sense to use a 
transmission range larger (in order sense) than the distance gained 
along the chain of home-points followed by the routing algorithm 
(otherwise the same transmission range could be used in a previous 
step to directly reach the destination, as done in the last step). 

The following lemma will be useful for dimensioning Ai . 

Lemma 2. At routing step i, the minimum value of Ai which 
guarantees that lirainin^ao Pi^ipopulated squarelet \ i} > is 
given by: 



Ai = i 



e(logn) 
Proof, (see Appendix lAt. D 



e 
e 



Zi 
2-i5 



<5< 1 

1<5<2 

5 = 2 
(5> 2 



(3) 



We observe that Ai depends both on the value Zi — 2^Zo of 
the associated routing step and on the power law exponent 5 which 
characterizes node mobility. 

For S < 2 the largest value of Ai corresponds to the last step 
(i.e., step 0, having minimum Zi equal to Zo). In this case, Zo 
cannot be chosen arbitrarily small. In particular, for any 5 < 1, 
condition \/!A7 = Ri = 0{Zi), coupled with ^, implies that 
Zo = Q,{n^'^). For 1 < 5 < 2, the same condition implies that 
Zo = 0, (n'^Y 

For S = 2, Ai does not depend on Zi, therefore the same trans- 
mission range R — ^/]ogn can be used in all steps. Notice that 
Ri — 0{Zi) is satisfied because we require that Zo — Q{\/\ogn). 

For S > 2 the first step (having the maximum Zi) is the one that 
requires the largest value of Ai. Also in this case, to apply Lemma 
[H we require that Zo — Q{^/\ogn). 

To summarize, the conditions on the free parameter Zo of our 



Zo = 



Q{n~ 



r2(Vlogn) 



5<1 

1 <5< 2 
5>2 



(4) 



At last, we need to dimension the probability distribution pf with 
which slots are assigned to the different steps of the routing algo- 
rithm. A natural choice is to equalize of the average number of 
transmissions that can occur for each step, and thus avoid that a 
particular step becomes the system bottleneck. This is obtained by 
making pf inversely proportional to Ni — -j-, which is an up- 
per bound to the average number of parallel transmissions E[Ni] at 
step i. 

Given that imax ~ Y\ log2 '^— log2 Zo\ is the maximum number 
of steps traversed by messages, using (|2} we have 



Nj 



A, 






(5) 



Ni 



We will later see that this choice ofpt indeed equalizes the average 
number of parallel transmissions occurring at each step. 

6.2 A first throughput characterization 

An upper-bound to the network throughput achievable by our 
scheduling-routing schemes can be easily computed in terms of the 
maximum number of messages that can flow in one slot from the 
sources to the destinations, under the optimistic hypothesis that 

Pr{active squarelet | i} = Pr{popuIated squarelet | i\ 

Considering a particular step i, the average number of messages 
stored at i-th relay nodes which can advance to (i— l)-th relay 
nodes (or be delivered to the destinations, in the case i = 0) in 
one slot is bounded by pfNi. It follows that an upper-bound to the 
network throughput A is expressed by: 



A = n \ = inm{piNi) 

i 

Given the design choice (O we obtain: 

A = poNo 



(6) 



J2i=o ^■<- 



For 5 < 1, being A, = 2_ll, it follows: 



A: 



For 1 < (5 < 2, being Ai = - 
A= ^ 



= e{n2Zo) 



(7) 



2-iS 



it follows: 



Q{n^Zo 



E^max n 2 

For S = 2, let ^0 = 71-^. We have: 



T "" 


/ Q(log^{n)) 


(3 < 1/2 


'T^imax 1 

T,^=o^ogn 


13 = 1/2 


it last in the case 5 > 2 
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(8) 
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6.3 Delay analysis 

Now we turn our attention to the delay. We consider a generic 
node a which is storing a message m in step i, directed to final 
destination d. We need to evaluate the average number of slots 
needed to make this message advance by one step. By so doing we 
neglect the effects of possible contention for transmission among 
different messages in step i stored at node a; this is equivalent to 
ignoring queueing delays within each node, and considering only 
the service time (we will analyze queueing effects in Section [6l4l l. 

Node a is enabled to transmit message m in a given slot if the 
following three conditions hold: i) the slot is assigned to the i-th 
step; ii) an eligible relay node b for the head-of-the-line message of 
the queue devoted to step i > 1, or a suitable destination node d 
(for step i = 0), is found in the same squarelet in which a resides; 
iii) among all the eligible contending pairs residing in the same 
squarelet, a pair is selected in which a acts as transmitter. The 
occurrence of condition i) is simply pf . The probability p" of the 
occurrence of condition ii) is computed in Appendix IB] 



to destination can be computed as Ds ~ X^j Dg. It follows: 



pT = < 



e(A^ 



e A,- 



2(1-5), 






5< 1 

1<5 <2 

S = 2 
5>2 



(11) 



In the above expressions, \Ti\ denotes the number of potential re- 
ceivers of message m that exist in the network. For i > 0, we 
have IFil = Q{Zf), whereas in the last step there exists a unique 
receiver (i.e., the destination), hence |ro| = 1. 

The probability p" of the occurrence of condition iii) depends on 
the value of S. For 5 < 2, pf = 0(1), because the average number 
^{^(ab)] of eligible contending pairs (a, 6) that can be found in a 
squarelet is finite (see Appendix IaI and by Jensen inequality we 
have pf > S[l/iVf,,)] > 1/S[<,)] = e(l). For S > 2, 
instead, there is an infinite number N^ of transmitting nodes a 
and a finite number N^ of receiving nodes b (see Appendix lAt. 
Therefore, under the pessimistic assumption that all nodes a have 
a packet to transmit, it follows pf > E[l/N^] > 1/E[N^]. We 
have, 

9(1) 
n I i°s" 

" \Ai log Ai 

n{i/Ai) ' 



pf 



(5< 2 
(5 = 2 
(5> 2 



(12) 



The probability p, that a transmits message m in a given slot can 
then be computed as pj — pfp"p^. 



Pi = i 



n2Z? 



e 



ir.l^S 



Z^ log^ (n) log log n 



(13) 



5< 1 

1<5<2 
5 = 2 
5 > 2 



At last, under the pessimistic assumption that all nodes in the net- 
work are constantly backlogged by messages in step i, the chances 
that message m is forwarded in a slot form a memoryless Bernoulli 
process, since pf depends only on the geometry of nodes, which 
completely regenerates at every slot. Thus, it immediately follows 
that the average service time D^^ of a message in step i is equal to 
l/pj ■ Neglecting queueing delays, the total delay Ds from source 



Ds 



Oi 



e{Zon2) S<1 

e{Z^n^-i) 1<5<2 

Zq log^(n) log log TI N A — 9 

log Zo ^ 



(14) 



0{n'^^-'Z^o) 



S>2 



6.4 The effect of traffic 

In the above computation we have assumed that the delay ex- 
perienced by a message at a node is of the same order of magni- 
tude as the service time, i.e., the time that elapses from when the 
message becomes head-of-the-line to when it is successfully trans- 
mitted. This can be justified by Kingman's upper bound to the total 
delay of a discrete-time GI/GI/1-FIFO queue [25], which states that 
as long as the second moment a^ of the number of simultaneous ar- 
rivals is finite, the second moment a"|j of service time is also finite 
and the queue-load p is strictly less than one, the average queueing 
delay D'' is bounded by: 



D' < D\ max 1 



ol 



2{D^sni^p) 

from which it follows that D' = O(Z)g). In our case the fact that 
(Ta < oo is immediate in light of the fact that at most one message 
can be enqueued per slot at any node, while cr^ < oo derives from 
the fact that the service time distribution is stochastically domi- 
nated by a geometric distribution. 

6.5 Maximum Throughput evaluation 

We observe that the whole system can be represented by an 
acyclic network of GI/GI/1-FIFO queues; indeed messages ad- 
vance in the network visiting queues associated to decreasing step 
indices i, which guarantees the absence of loops (see Figure |6). 
As a consequence, offered traffic A can be successfully transferred 
through the network, as long as no queue is overloaded. 




Figure 6: Internal architecture of tlie nodes: queue structure 

The stability condition for the single FIFO queue present at each 
node a, which stores messages belonging to step i > 1 of the rout- 
ing algorithm, is Af < pj , where A" is the aggregate arrival rate 
of messages in step i at node a, while pf is the transmission prob- 
ability computed in Section |63] 



Considering the last step (i=0), since in this case messages di- 
rected to different destinations are stored in separate queues, the 
stability condition for the generic queue storing messages directed 
to destination d is Ao"* < po , being Af}'' the arrival rate at node a 
of last-step messages directed to d. 

Substituting the expression for pj , we immediately obtain that 
the maximum traffic sustainable by node queues satisfies: 



xt = e(A) 



i > 1 



and 



\ ad 
^0 



e 



i.e., every node a can sustain, for each step, a traffic B(A), ex- 
pressed in messages per slot (in the case of step i = 0, recall that 
every node w.h.p. delivers last-step messages to Q{Zq) different 
destinations). 

Note that our scheduling-routing scheme w.h.p. distributes the 
network traffic uniformly (in order sense) among all network nodes. 
We conclude that the sustainable network throughput is w.h.p. 
0(A), i.e., the upper bounds previously computed in Section [6^ 
are asymptotically tight. 

We remark that this last finding implies that: i) no step becomes 
the system bottleneck under our choice of slot assignment probabil- 
ities pf; ii) if the network sustains a throughput 0(A), necessarily 
Pr{active squarelet | i} — B(Pr{populated squarelet | i}). 

At last, we can claim the optimality of our scheduling-routing 
schemes for 5 < 2. This in light of the following proposition. 

Proposition 1. Consider any scheduling-routing scheme ac- 
cording to which messages are delivered to destinations d by nodes 
a whose home-points satisfy d^^ = Q{Z), employing a trans- 
mission range R. Then, necessarily the network throughput is 
A = 0{-^). In addition, if no message replication is employed, 
the delay is D = Q. [1/Pr{meet}), where Pr{meef} is the proba- 
bility that any two nodes (a, d), with d^d = 0(-^). fall within the 
transmission range of each other. 

Proof. The throughput condition is an immediate consequence 
of the protocol interference model, that prevents two nodes whose 
relative distance is smaller than -^ R from transmitting simultane- 
ously [1]. The delay condition immediately descends from the ob- 
servation that a can deliver a message to d, only if d falls within 
the transmission range of a. D 

Comparing the bounds of Proposition [T] with the performance 
achievable by our scheduling-routing schemes (for which Z = Zo 
and R — ^/Aq), we conclude that our schemes, for 5 < 2, achieve 
optimal delay-throughput trade-offs within the class of algorithms 
that do not employ message replication. 

6.6 Possible delay throughput trade-offs 

Let us first examine the case S > 2. According to JlOt the sys- 
tem throughput does not depend on Zo, whereas the delay l ll4t in- 
creases with Zq. Thus, Zo should be minimized to reduce as much 
as possible the delay, while no delay-throughput trade-offs are pos- 
sible. Recalling the condition Q we set Zo — ^/logn and obtain 



X 



Also for 5 — 2 the dependence of the throughput on Zo is weak, 
since any choice of Zo = 0{nP), with /3 < 1/2, leads to the same 
throughput, while a log n factor is gained selecting Zo ~ 0(n^ "), 
however at the expense of a severe increase of the delay. Hence it 
appears optimal to select Zo ~ 0\/logn, resulting in: 



A = e 



log-'(n) 



D = 0{log''{n)) 



For S < 2, instead, different trade-offs between throughput (opti- 
mized selecting Zo large) and delay (minimized selecting Zo small) 
are possible. Let Zo = n^. 

For 1 < 5 < 2, being |f^ < P < 1/2, we obtain 

A = e(n*/^-^+«2-^)) D = e(n^-*(^/^-'^') 

whereas for 5 < 1, being 1/6 < /3 < 1/2, we obtain 

X = e{n^-i) D = e(n'^+3) 

6.7 Alternative scheme for large s 

It should be noticed that the performance of our scheduling- 
routing schemes in the case of S > 2 rapidly deteriorates in terms 
of both throughput and delay for large values of S. Therefore, we 
propose here an alternative scheme, valid for 5 > 2, whose perfor- 
mance does not depend on 5. 

We observe that, for 5 > 2, every node spends a constant fraction 
of time within a finite distance from its home-point. Hence it is 
possible to devise a scheduling-routing scheme which is based on 
this property and that does not exploit node mobility at all. The 
network area O is divided into a regular square tessellation whose 
elements have area equal to log n. Within each squarelet, only the 
nodes whose home-points belong to the squarelet itself, are allowed 
to transmit or receive, using a fixed transmission range equal to 
yTogri. Since each node spends a constant fraction of time in the 
squarelet containing its home-point, the performance that we can 
achieve is the same as if nodes were fixed at their home-points, i.e., 
it is equivalent in order sense to the Gupta-Kumar case. We can 
therefore obtain a per node- throughput A = 0{l/\/nlogn) and 
delayi3s = e(Vw/logn)[8]. 

A comparison with our previous scheduling-routing schemes 
leads to the conclusion that the alternative scheme becomes con- 
venient for 5 > 3 (see Figure[T]l. 

7. THE SLOW MOBILITY CASE 

In presence of slow mobility it is potentially possible to deliver 
messages along multi-hop paths in a single slot. Our goal is to 
understand whether this possibility can be exploited within our 
class of scheduling-routing schemes to further improve their per- 
formance. 

A particular case, which serves as a lower bound on the sys- 
tem performance, is when there is a single step, characterized by 
Zq — v'n, for which the resulting scheme essentially degen- 
erates to the Gupta-Kumar case, providing per-node throughput 
A = f2(l/\/nlogn) and delay Ds ~ 0(1) (according to the new 
definition of slot). 

Potentially better delay-throughput trade-offs could be achieved 
by employing a hybrid scheme in which messages are first routed 
according to the previously described bisection routing scheme, up 
to arriving at a critical distance from the destination (to be spec- 
ified), and then transferred to the destination in just one slot over 
a multi-hop path. This means that when the last step (i = 0) is 
scheduled, the network area is subdivided into squarelets of area 
Bo ~ uj{Ao), and many (disjoint) multi-hop transmissions are es- 
tablished within each squarelet. This class of hybrid schemes may 
be effective for S < 2., i.e., when performance is indeed limited by 
the last step. Therefore in the following analysis we assume S < 2. 

7.1 Analysis 

We first analyze the throughput performance in the last step, con- 
sidering the effect of multi-hop transmissions. Denoting by Q{Bo) 
the number of multi-hop transmissions that can be enabled in a 
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Figure 7: Power A/D (in log scale) as a function of S (slow mo- 
bility case). 



squarelet of size Bo during a slot, and by A'o the total number of 
parallel transmissions occurring in tire network at step 0, we liave 



E[No] = ^E[Q{Bo)] 



(15) 



In the previous bisection scheme, the squarelet size Ao was dimen- 
sioned in such a way that a finite number of eligible transmitter- 
receiver pairs (single-hop) can be found in a squarelet of area Aq. 
Therefore, in a squarelet of area Bo we have Co = Q{Bo/Ao) 
nodes that can communicate among them over multi-hop paths. 
Although the number of pairs that can potentially be established 
increases quadratically with Bo, only Co of them are not conflict- 
ing (i.e., do not have a node in common). Therefore, we have 
E[Q{Bo)] = Co = Q{Bo/Ao). Plugging this expression in ^ 
we immediately see that E[No] = -E[A^o] = n/Ao for any choice 
of Bo, hence there is no advantage in terms of spatial reuse by using 
multi-hop paths. 

Turning our attention to the delay of the last step, we consider a 
message m stored at a and destined to d. The probability po that 
the pair of nodes (a, d) reside in the same squarelet of area Bo is 
(see appendix IB}: 



Po 



e(^) 

72(l-a) 

e(Bo%^ 



5< 1 
1 <5 <2 



hence it increases linearly with Bo. However in this case the proba- 
bility pp that (a,d) is selected is O(l/Co) = 0(Ao/Bo), because 
among all 0(Co ) pairs that can potentially be enabled only 0(Co) 
are not conflicting. As a result, the product poPo does not depend 
on Bo, therefore no gain in terms of delay is obtained by a scheme 
that attempts multi-hop transmissions at the last step. The only ex- 
ception to previous arguments is represented by the case in which 
there is a single step, and thus the hybrid scheme degenerates to the 
Gupta-Kumar case. 

A graphical representation of the network power achievable in 
the slow mobility case is reported in Figure |7] We observe that 
our class of routing-scheduling schemes performs better than the 
degenerate case of Gupta-Kumar for values of 5 comprised in the 
interval [1.6,2.5]. 

8. DISCUSSION 

Although mainly of theoretical interest, we believe that our work 
can provide fundamental principles to design smart routing proto- 
cols for mobile ad-hoc networks exploiting the spatial diversity of 
the nodes. In particular, our bisection routing strategy appears to be 
an attractive solution as long as the spatial distribution of users fol- 
lows a power law. Although there is strong experimental evidence 



that this is true in many context related to both human and vehicular 
mobility, the precise exponent of the power law needs to be care- 
fully estimated, either through measurements or by an autonomous 
self-learning procedure. Indeed, two main regimes occur depend- 
ing on the exponent value: for values smaller than 2, the larger de- 
lays are expected in the last hop, whereas above 2 the critical hop 
is the first one. To avoid large delays, this requires to increase the 
transmission range either at the end or at the beginning of the route, 
which is a rather new concept. In this paper we have not addressed a 
number of issues which are essential to translate our scheme into an 
implementable solution. In particular, the impact of more realistic 
mobility models, and the local dissemination of information about 
the home-points of the nodes. However, we believe that the above 
issues do not compromise the applicability of our results to a real 
setting. Finally, the idea of progressively narrowing the selection of 
relay nodes through a divide-and-conquer technique, based on the 
similarity of the mobility pattern of the nodes, could be an interest- 
ing design principle to apply in a wider sense than presented here. 
In this case, one has to evaluate the mixing degree of the nodes 
in terms of distance scale membership, having defined a suitable 
distance metric between the mobility pattern of two nodes. 

9. CONCLUSIONS 

Previous work on delay-throughput trade-offs in mobile ad hoc 
networks has mostly been done under the assumption that nodes 
are homogeneous and uniformly visit the network area. In this pa- 
per we have shown that this condition can be largely suboptimal. 
Restricted mobility, which is usually found in practice, can be ex- 
ploited by intelligent scheduling-routing schemes which make use 
of the geographical information about the location most visited by 
a node. In particular, we have introduced a new class of schedul- 
ing routing schemes which significantly outperform all previously 
proposed schemes for a wide range of restricted mobility patterns. 
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of pairs averaged over all possible instances of H. Hence, we have: 



2 



(t,{Ha)4>{Hi)I{Ha, Hb) AHa dHb = 



Al 



(j>{Ha)<f}{Hb)I{Ha, Hb) AHa AHb 

Haeo ^H"[,eci:||ir6||>||Hall 



= Al / ct>{Ha)<t>{Hb)I{Ha,Hb)AHaAHb + 

+ Alf f (j,{Ha)<f){Hb)I{Ha,Hb)AHaAHb 

where the last two terms, for large C, may be interpreted, respec- 
tively, as the contribution of pairs (a, b) whose home-points are at 
distance 6(Zo) and the contribution of pairs whose home-points 
are at distance uj{Zi3). 

The first term can be approximated by the integral 



J H, 



AU I <j)\Zo)I{Ha,Hb)AHaAHb 



Q[Al<j,\Zo)Zt) 



since, by triangular inequality, it must hold that: 

lli^bll < \\Ha-Hb\\-\\Ha\\={C+l)Zo 

The second term instead provides a contribution 



Q[Alzl 



Hc,&0:\\Hcu\\>CZo 



\Ha)AHa) 



APPENDIX 

A. DIMENSIONING OF Ai 

We first focus on the last step (i = 0). We consider all distinct 
pairs of nodes (a, h) whose home-points satisfy Zo/2 < \\Ha — 
Hb\\ <iZo/4.. 

Given a squarelet S of area Ao, whose center of mass is for sim- 
plicity placed at the origin, we first evaluate the average number 
E[Nf^ i))|H] of 'far' node pairs (a, 6) (i.e., pairs such that both 
||-ffa|| > ,^o/4 and IIHftll > Zo/4) that fall within S. We further 
assume \/Ai < cZq for some c < 1/4. 

Let us denote with I{Ha,Hb) the function that returns 1 if 
Zo/2 < \\Ha - Hb\\ < 3Zo/4 and both \\Ha\\ > Zo/4 and 
\\Hb\\ > Zq/A. By construction, it follows: 



^ a bjta-'XiesJx2es 

(l>{Hb - X2)I{Ha, Hb) dXi dXa 



Substituting the expression for <j>{X), we obtain that the number of 
'far' node pairs is; 



E[Nfa,b)] 



^ n ■' 
.2^2(2-5) 

"■ .9 

e( 



log^ n ' 



I Q{AiZt, 



2r,2Ci-S), 



&< 1 

1<5<2 
5 = 2 
5 > 2 



To ensure that a far pair (a, 6) is found within S with a non- 
vanishing probability it is necessary to dimension ^o in such a way 
that E\Nf^ (,)] is not vanishing as n ^ oo. On the other hand mak- 
ing E\Nf^ 5)] = 0(1) is also sufficient to guarantee that a far pair 
is found within squarelet S with non vanishing probability. 

Now we turn our attention to 'near' pairs, i.e., pairs of nodes 
(a, &) in which one of the two nodes (let it be node a) has home- 
point satisfying: \\iia\\ < Zq/A. First, note that by triangular 
inequality it must be ^o/ 4 < \\Hb\\ < Zq .Thus to find a 'near' pair 
(a, b) within 5 with finite probability, we must find within 5 at 
least one node a whose home-point satisfies \\Ha\\ < Zq/A and at 
least one node b whose home-point satisfies Zq/A < \\Hb\\ < Zq. 

The average number E[Na] of near nodes a in squarelet 5 can 
be evaluated as; 



Note that E[Nf^ f^AH] is a random variable over the space 
of all the possible home-point locations H. However, since 
Zo — r2(yiogn), applying Proposition [T] it can be immediately 
proved that, for n ^ oo, with high probability, E[Nf^ 6)|H] = 
e(£[iVf,,,)]), being i?[iV(i,,)] = E»[E[Nf^^,^] \ H] the number 



S[iVa]xA 



/ (l>{Ha) AHa 

JHaeO:\\Ha\\<^ 



e(. 



■^0^0 



e(Ao) 



S <2 



©(^oi^) 5 = 2 



S>2 



whereas the average number E[Nb] of nodes b within 5, such that 



l/4Zo < \\Ht\\ <5/4Zo,is 



E[Nt 



J H, 






n 2 
^ log n ' 

l Q{AoZl-' 



-) 5<2 

5 = 2 
5>2 



In order to have near pairs within squarelet S with a non vanishing 
probability, both E[Na\ and E[Nh\ must be non vanishing. At last, 
Ao must be dimensioned in such a way that either both E[Na\ and 
E[Nb\ are non-vanishing, or E[Nf^ ,,>] is non-vanishing, leading to 
the expressions in (O. for i = 0. 

Now turning our attention to the case i > 0, we have to evaluate 
the average number of distinct pair of nodes (a, 6), whose home- 
point satisfy |^i < 1 1 //a — ^i, 1 1 < ^Zi, falling within a squarelet 
S of area Ai, whose center of mass is for simplicity placed at the 
origin. 

This can be done using exactly the same approach as before, i.e., 
distinguishing between 'far' node pairs (a, b) (pairs such that both 
\\EIa\\ > Zi/IQ and \\Hb\\ > Zi/16) and 'near' node pairs (pairs 
(a, b) in which at least one of the two nodes (let it be node a) has the 
home-point satisfying \\Ha\\ < Zi/16), and evaluating separately 
the two contributions. 

B. EVALUATION OF Pf 

Consider a pair of nodes a and b such that 1 1 Ha — Hb \\ = D. We 
wish to evaluate (in order sense) the probability P(a,b) that (a, 6) 
are found within the same squarelet S of area A: 



P(a,b) 



^Pr{aG5}Pr{6G5} 
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4>{X2 - Hb) dXa dXi 



X2&S(Xi) 



4,{Xl-Ha) 0(X2 

XiGO Jx2SS(Xi) 



Hb) dXa dXi 



being S{Xi) the squarelet containing Xi. Let V C O denote the 
set of points X closer to //a than //ft (i.e., such that: ||i/a— X|| < 
\\Hb — X||), then by symmetry: 



P(a,b) 

= 2 



2 / (l>{Xi - Ha) / 0(^2 - Hb) dXa dXi 
iXiev Jx2es{Xi) 



(j}(Xi - Ha) / <t>(X2 - Hb) dXi dXi + 
Xiev»:|lXi-H„||<£ Jx2eS(,Xi) 



+ 2 



XiSV" 



0(Xl - Ha) 
-<nXi-H„||<4D 



(P(X2-Hb)dX2dXi+ 

X2es(Xi) 



+ 2 / ?i(Xi - Ha) / 0(^2 - Hb) dXa dXi 

J XieV'^:\\Xi-Ha\\>4,D J X2eS{Xi) 

Now assuming VA < c D, being c <C 1, the first term in the above 
expression provides a contribution 



q(a<J}{D) j ^</.(Xi)dXi) 



XilK 



since by triangular inequality on the considered domain, at the same 



Ha-XrW >\\Ha-Hb\\- 

Hb-Xi\\ < \\Ha-Hb\\ + 
X2-X1II <cD. 



time||//i,-Xi||> WHa-HbW- 
i/a||-||Xi|| > 3D/4-cDand 
Ha-Xi\\ < 5D/4 + cD, while 
The second term provides instead a contribution 
O (ylD^0^(Z))) since on the considered domain it must be 
7:^/4 < \\Hb-Xi\\ <5D + c. 

At last the third term provides a contribution 



e A 



J XieV<':\\Xi-Ha 



\\>4D 
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since on the domain \\Ha — Xi\\ = Q{\\Hb — X2I 
Substituting the expression for (i>(X) we obtain: 



P{a,i,) = < 



e(-) 

r,2(l-5x 



<5<1 

1<<5<2 

5 = 2 
(S> 2 



At last, observe that 'P(a.b) equals po if we set 
Zo/2 < D < 3^0/4, and either A = Ao in case of fast 
mobility or A = _Bo in case of slow mobility. 

Now, to evaluate the expression of pf for i > 1, let F, denote 
the set of all potential nodes b which can receive a message in step 
i from node a. By definition: 



P? ^Yl P^^'' ^ ^^ P^{^^ G r. : fe G 5} = 

= J2^r{aeS}[l-ll{l-Pr{beS}) 

s beFi 

In addition, since it can be proved (as shown below) that: 
1 - n (1 - P^i'^ e "5})) = e(min(l, J2 Mfe G 5})^ 

(16) 



beVi 
we can conclude that p° 



= e I V / <^(Xi - Ha) mm( I, V / <7i(X2 - Hb)dX2)dxA 
\ s ixies ^ ber.-'^^es ' j 

Then operating in a way similar to the computation of P(a.b) we 
obtain the expressions in dllb . 

To conclude the proof we have to show how to obtain \\b\ . The 
proof if trivial if there is at least one b* G Vi such that Pr{fe* G 
5}) = 0(1), since by construction 



1 - n (1 ^ P^'^^ ^ ^^n > P^'^^* ^ "^^ 



So, suppose that Pr{6 G 5}) = o(l) for all 6 G F^. In this case, 
first note that 

log n (1 - Pr{& e 'S}) = E log(l - P^{^ '^ S^') 

bsTi ber; 

Then, since —2a; < log(l — a;) < x for any < x < xq, with 
Xf) ~ 0.8 solution of equation 1 — xo = exp(— 2xo), it follows: 

-2^ Pr{6 G 5} < log f]^ (1-Pr{6 G 5}) < -^ Pr{6 G 5} 

ber; ber, ber; 

for X^bgr Pr{fe G 5} < xo, from which: 

1 - exp(- Y^ Pr{6 G 5}) < 1 - H (^ ~ ^^^^ ^ '^}) ^ 

< 1 - exp(-2 Y Pi'i^ G '^}) (1'^) 

ber; 

Now 1 - exp(- Eber, Pi^{& G 5}) > 1/2 Eter. Pr{b G 5} 
for Ebgr P^{^ G 5} < Xo, thus 

1/2 Y P^{^ G 5} < 1- J]^ (1-Pr{6 G 5}) < 2^ Pr{& G 5} 

ber,; ber; ber, 

i.e., the assertion is proved for Eber Pr{6 G 5} < xq. The 
extension to the case Etgr Pr{b G <S} > xq is trivial in light 
of the fact that (1 — Hber (1 ~ Pr{& G 5})) is increasing with 
respect to Pr{fe G 5}. 



